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SECTION  I 


INTRODUCTION 

This  report  contains  tables  of  probabilities  of  damage  to 
point  or  circular  normal  area  targets  using  the  cumulative  log- 
normal damage  function.  The  probability  of  damage  for  such  a target 
is  a combination  of  three  probabilities:  the  probability  that  a target 

receives  at  least  a certain  amount  of  damage  (called  the  distance 
damage  function);  the  distribution  of  weapon  impact  points  (CEP);  and  the 
density  of  the  target  area  (R95).  In  order  to  easily  evaluate  the  com- 
bined probabilities,  we  utilize  analytic  approximations  for  these  sep- 
arate distributions.  The  complement  of  the  cumulative  lognormal  func- 
tion is  used  as  tiie  analytic  approximation  to  the  distance  damage 
function,  and  the  circular  normal  function  represents  the  joint  density 
of  the  impact  distribution  and  the  target  density. 

Prior  to  1972,  the  distance  damage  function  was  described  by 

a circular  normal  distribution.  The  joint  distribution  of  the  distance 

damage  function,  impact  distribution  and  target  density  was  referred  to 

as  the  circular  coverage  function.^  Tabular  evaluations  of  the 

circular  coverage  function  had  been  in  use  for  many  years.  However, 

since  1972,  the  complement  of  the  cumulative  lognormal  function  has 

3 

been  used  to  approximate  the  damage  function.  In  adapting  this  approx- 
imation for  use  at  Science  Applications,  Inc.,  adequate  evaluations  for 
comparison  and  verification  could  not  be  located.  Consequently,  these 
tables  were  produced  for  this  purpose.  The  following  pages  give  a brief 
description  of  the  geometrical  meaning  and  mathematical  background  of 
the  tables  and  a brief  discussion  concerning  the  formulation  and  use  of 
the  tables.  A more  detailed  explanation  of  the  mathematical  background 
and  programming  techniques  used  to  generate  these  tables  can  be  found 
in  Reference  3. 

* ' H . H . Germond,  "The  Circular  Coverage  Function",  Project  RANI) 

Research  Memorandum  RM-330,  The  Rand  Corporation,  26  January  1950. 

2 "Offset  Circle  Probabilities",  R-234,  The  Rand  Corporation, 

14  March  1952. 

3- "Mathematical  Background  and  Programming  Aids  for  the  Physical 
Vulnerability  System  for  Nuclear  Weapons",  DI -550-27-74 , 

Defense  Intelligence  Agency,  1 November  1974. 
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SECTION  II 


PROBABILITY  OF  DAMAGE  ANALYSIS 

This  Section  gives  a geometrical  interpretation  and  a mathe- 
matical definition  of  the  problem  at  hand.  It  is  drawn  largely  from 
Reference  3 and  unpublished  material  by  the  authors. 

The  basic  problem  is  to  calculate  the  probability  of  tar- 
get damage  from  a single  weapon.  If  the  target  has  an  extended 
area,  the  probability  of  damage  is  replaced  by  the  expected  fraction 
of  the  target  damaged.  In  this  report,  the  term  "probability  of  damage" 
is  used  for  both  point  and  extended  area  targets.  The  solution  of  this 
problem  requires  the  integration  of  three  independent  probability  dis- 
tributions— the  distance  damage  function,  the  weapon  impact  distribu- 
tion and  the  target  density  distribution.  The  distance  damage  function 
is  represented  by  the  complement  of  the  cumulative  lognormal  function. 
The  weapon  impact  distribution  and  the  target  density  are  represented 
by  the  circular  coverage  function.  A geometric  interpretation  follows: 

Assume  that  a weapon  is  aimed  at  point  A,  which  is  offset  from 
the  target  T the  distance  D,  as  shown  in  Figure  1.  The  distance 
damage  function  is  defined  by  the  weapon  radius  (WR)  and  the  distance 
damage  sigma  (o^)  (see  Figure  1A).  WR  is  the  radius  of  a circle 
centered  at  the  impact  point  that  contains  as  much  area  inside  that 
is  undamaged  to  a specified  level  as  there  is  damaged  outside.  a ^ is 
the  standard  deviation  of  the  distance  damage  function  divided  by  the 
weapon  radius.  Analytically,  we  approximate  P^,  the  distance  damage 
function  by  the  complement  of  the  cumulative  lognormal  function  as 
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Figure  1. 
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Using  the  distance  damage  function,  we  can  calculate  an  average 
probability  of  damage  to  a point  target  for  each  detonation  point  by 
weighting  P()  by  the  probability  that  the  weapon  lands  at  that  detona- 
tion point,  called  the  impact  distribution. 


The  circular  normal  impact  distribution  is  depicted  graphically 
in  Figure  IB.  It  is  defined  by  the  Circular  Error  Probable  (CEP),  which 
is  a circle  about  the  aimpoint  which  contains  50  percent  of  the 
impacts.  The  probability  that  a weapon  will  impact  within  a distance 
R<  r <^R  + AR  of  a target  T , assuming  an  aimpoint  offset  of  0 is: 


f(r) 


2tt 

f rV 

J0  2naCEP 


-(p2/2°CEp) 


de  R AR 


R2  + l)2  - 2 HI)  cos  -) 


° CEP  = CEP  / In  2 


If  we  have  an  extended  area  target,  we  must  consider  the  target 

density  distribution.  It  is  illustrated  in  Figure  1C  and  is  also 

assumed  to  be  a circular  normal  distribution.  It  is  defined  by  the 

R95  radius  which  is  a circle  centered  at  the  target  which  contains 

95  percent  of  the  target  value  distribution.  Since  the  target 

density  f(t)  and  the  delivery  error  f(r)  are  independent  normal 

distributions,  the  joint  distribution  is  f(r,t)  = f(r)  • f(t),  and  it 

is  referred  to  as  a circular  normal  distribution.  Also,  the  variances 

2 2 2 2 

can  be  combined  as  a - a.  + a„T,T,,  where  a is  the  variance  of  the 

a t CEP  a 

joint  distribution  . Hence, 

? 1 n 2 ? 1/^ 

CEPa  = (CEP2  +TlrT5o-  R95  } • 


To  find  the  probability  of  achieving  the  desired  level  of 
damage,  PD  , to  the  specified  target  T , we  must  integrate  all 
three  distributions: 

PD  = 


i 


pd(r) 


f ( r ) dr  . 
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SECTION  III 


TABULAR  RESULTS 

This  section  explains  how  to  use  the  tables.  It  includes  a 
nomograph  which  tells  when  and  how  to  replace  CEP  with  the  adjusted 
CEP  (CEP, ) in  order  to  account  for  distributed  area  targets.  It  also 
explains  the  organization  of  the  tables  in  terms  of  the  parameters  used 
and  gives  several  examples.  Finally,  it  discusses  the  table  limits 
as  certain  variables  increase  without  limit. 

3.1  USE  OF  ADJUSTED  CEP 

No  CEP  adjustment  is  necessary  for  point  targets.  For  cir- 
cular normal  distributed  area  targets,  the  target  density  distribu- 
tion can  be  included  by  adjusting  the  CEP  as  explained  in 
Section  2.  The  exact  equation  for  the  adjusted  CEP  is: 

cl 

f 9 I n ° 9 1/2 

CEP ( = CEP  + yip ( R95  * 6076.1155) 

r 2 2i1/'2 

= I CEP  + 8542294  * R95  J 

The  "In  2/ In  20"  term  converts  the  95th  percentile  of  the  R95  dis- 
tribution to  the  50th  percentile  used  for  CEP.  The  "607(5.1  155"  factor 
converts  nautical  miles  to  feet. 

As  can  be  seen  by  the  equation,  CEP  is  composed  of  two  terms, 

CEP  and  R95.  Usually  one  or  the  other  of  these  terms  dominate. 

Figure  2A  shows  the  limits  of  CEP,  for  a wide  range  of  values  of  CEP 

and  R95.  The  diagonal  line  is  for  an  R95  of  zero  and  CEP  equals 

CEP.  The  vertical  lines  are'  for  when  R95  dominates  and  CEP  equals 

R95  imes  2923.  Also  indicated  are  nomographs  (Figures  3 and  4)  from 

which  one  can  obtain  more  precise  values  of  CEP,  , given  CEP  and 

u 

R95.  Figure  213  shows  when  to  use-  the  nomographs  and  when  to  use  the 
approximations.  The  boundaries  were  chosen  to  ensure  an  accuracy  of 
at  least  .01  in  the  final  probability  of  damage  (PD)  result. 
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3.2 


TABLE  ORGAN  I ZAT ION 


As  explained  in  Section  2,  in  order  to  determine  the  proba- 
bility of  damage  (PD)  to  a pointer  circular  normal  area  target,  it 
is  necessary  to  integrate  over  three  distributions.  These  distribu- 
tions are  represented  by  the  following  parameters: 

Distance  Damage  Function:  WR, 

Weapon  Impact  Distribution:  CEP,  D 

Target  Density:  R95 

As  mentioned  above,  the  CEP  and  R95  can  be  combined  as  the  adjusted 
CEP&  , leaving  us  with  four  parameters  (WR,  , CEP  , D)  . By 
normalizing  one  of  the  parameters  and  by  using  separate  tables  for  dis- 
crete values  of  we?  reduce  the  presentation  format  to  double-entry 

tables  of  PD  . 

3.2.1  Standard  Parameter  Set 

The  probability  of  damage  (PD)  to  a point  or  circular  normal 
area  target  has  been  evaluated  as  a function  of  WR/CEP,  versus  D/WR  . 
WR  is  the  weapon  radius,  CEP^  the  circular  error  probable  of  the 
joint  dist r ibut ion  of  the  target  density  and  the  impact  distribution, 
x WR  the  standard  deviation  of  the  distance  damage  function,  and 
D the  offset  distance  from  the  aimpoint  to  the  target.  The  standard 
data  set  is  chosen  as  follows: 

a.:  10,  20,  30,  40,  50%  of  WR 

CEP  : 1000  feet 

a 

WR/CEP  : 0.0  - 20.0  in  varying  intervals 

(I 

I)/WR:  0.0  - 4.8  in  varying  intervals. 

PD  has  been  tabulated  in  terms  of  WR/CEPo  to  eliminate  CEP  as  a 
table  input  variable.  The  intervals  have  been  chosen  so  that  linear 
interpolation  between  tabulated  values  will  give  three  digit  accuracy 
(±.002).  Interpolation  between  tables  of  o.  is  not  recommended  as  the 
accuracy  is  unknown.  The  upper  bounds  on  WR/CEP, ( and  D/WR  were  chosen 
to  exhibit  PD  approaching  definite  limits.  These  limits  are  dis- 
cussed in  Section  3.4. 
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3.3  USE  OF  TABLES  - EXAMPLE  PROBLEMS 

3.3.1  Example  1 - Point  Target 

A single  bomb  with  a weapon  radius  of  2500  feet  is  aimed  at  a 
point  2000  feet  from  a point  target . The  weapon  has  a 1000  foot  CEP 
accuracy  and  the  standard  deviation  of  the  distance  damage  function 
sigma  is  thirty  percent  (30%)  of  the  weapon  radius.  Hence, 

CEP  = 1000  feet 
D = 2000  feet 
WR  = 2500  feet 
o d = 30%  of  WR 
R95  = 0 

CEP^  = CEP  = 1000  feet 
WR/CEPa  = 2500/1000  =2.5 
D/WR  = 2000/2500  =0.8 

becking  the  table  for  = 30%  of  WR  (page  36),  we  obtain. 

PD  = .5572 

3.3.2  Example  2 - Area  Target 

A weapon  is  aimed  at  a point  8000  feet  away  from  the  center  of 
a circular  normal  distributed  area  whose  R95  radius  is  one  nautical 
mile.  The  weapon  radius  is  6.4  Kft  and  the  distance  damage  function 
sigma  is  forty  percent  (40%)  of  the  weapon  radius.  The  CEP  is  1300  ft 
Hence  , 


CEP  = 1300  feet 

D = 8000  feet 

WR  = 6400  feet 

od  = 40%  of  WR 

R95  = 1 n . mi . 

CEP  = 
a 

T 2 2 1 

[1300  + 8542294(1  )J  = 3200  feet 

WR/CEPa  = 

6400/3200  =2.0 

D/WR  = 8000/6400  = 1.25' 
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Checking  the  table  for  = 40%  of  WR  (page  46),  we  obtain: 

PD  = .2201 

3.3.3  Example  3 - Table  Interpolation 

In  example  1,  assume  that  the  weapon  radius  is  changed  to 
3000  feet  and  the  CEP  to  900  feet.  Hence, 


CEP 

a 

WR/CEP^ 

D/WR 


CEP  = 900 

feet 

D = 2000 

feet 

WR  = 3000 

feet 

od  = 30% 

of  WR 

R95  = 0 

CEP  = 900 

feet 

3000/900  = 3 1/3 
2000/3000  = 2/3 


To  find  PD  we  must  use  double  linear  interpolation  between 
the  table  entries  from  page  38. 


WR/CEP 

(1 

3.20 

3 1/3 

3.40 

C1 

B 

C2 

. 66 

. 7214 

. 7320 

D1 

T11 

B1 

T12 

D/WR 

2/3 

A 

A1 

PD 

A2 

.68 

. 7037 

. 7139 

D2 

T21 

B2 

T 

22 

14 


Given  the  data  in  the  table  above,  there  are  two  equivalent 
ways  to  find  the  desired  answer,  PD.  We  can  solve  for  A^  and  A,, 
and  then  interpolate  for  PD  , or  we  can  solve  for  and  B^  and 

then  PD.  Both  methods  are  shown  below. 


Interpolation  Method  A 

A1  ' T1 1 * <T21-TU>  [<A-D1)/(D2-D1j] 

= .7214  + ( .7037-. 7214)[. 00667/. 0200] 
= .715 

A2  - T12  * <T22-T,2)[<A-D1)/(D2-d1)] 

= .7320  + ( .7139-. 7320) [.00667/. 0200] 
= .726 

PD  = Aj  + (A2-A1)  [(B-C1)/(C2-C1)] 

= .715  + (. 726-. 715) [. 133/ .200] 

= . 722 


I nterpol at  ion  Method  B 

B1  ■ T1 1 * <T12“T11 )/<C2"Cl J] 

= .7214  + (. 7320-. 7214) [. 133/ . 200] 

= . 728 

B2  ■ T2]  * (T22-T21>[<B-C1>/<C2-C1>] 

= .7037  + ( .7139-. 7037) [. 133/. 200] 

= . 710 

PI)  = Br  + (B2-B1)f(A-D1)/(D2-D1)] 

= .728  + (. 710- . 728) [. 00667/ . 0200] 
= . 722 
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3 . 3 . ‘1 


Example  4 


- Required  Offset  Distance  for  Given  PD 


In  example  2,  assume  that  it  is  desired  to  restrict  the  damage 
to  the  area  target  to  no  more  than  ten  percent  (10%).  What  is  the 
required  aimpoint  offset  distance? 


CEP 

a 

WR/CEP^ 


CEP  = 1300  feet 

WR  = 6400  feet 

a , = 40%  of  WR 
d 

R95  = 1 n . mi . 

PD  < 0.10 

T 2 2 i1/2 

1 1300  + 8542294(1  )J  = 3200  feet 

6400/3200  = 2.0 


We  must  now  interpolate  in  the  table  on  page  46  to  find  the  value 
of  D/WR  which  will  give  PD  = 0.10. 


PD 

. 1120  „ 

. 1000  _ 

.0991  T,  1 

pi 

A 

P*l 

D/WR 

1.55  ^ 

99 

da 

1-60  d2  ! 
’ 

DA  = 

Di  + <W 

[(Pfl-Pp/(p2-P 

*>] 

= 1.55  + .05  [( .10-. 1120)/( .0991-. 1120)] 
= 1.5965 


D = D.  * WR  = 1.5965  * 6400  = 10218  feet 
A 
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3.4 


LIMITS  OF  TABULAR  EVALUATIONS 


As  mentioned  in  Section  3.2.1,  the  upper  hounds  on  the 
table  entries  for  WR/CEP  and  D/WR  were  chosen  to  show  PD  approach- 
ing  definite  limits  as  WR,  CEP,  , and  D inrrea.se/deerease  without 
limit.  In  particular,  the  limit  of  PD  as  WR/CEP  goes  to  iniinity 
is  the  distance  damage  function  and  the  limit  of  PI>  as  D/WR  goes 
to  infinity  is  zero.  A geometrical  explanation  ol  these  limits  is 
given  below. 

3.4.1  Limit  of  PD  Along  Horizontal  Axis 

Consider  the  variation  in  PD  as  WR/CEP^  increases  while 
D/WR  remains  constant.  Assume  CEPi(  is  a constant  , so  only  WR 
is  increasing  as  WR/CEP  increases.  Since  D/WR  maintains  a constant 

cl 

value,  D is  also  increasing  as  WR  increases.  Since  D and  WR  are 
both  increasing  while  CEP  remains  constant,  eventually  the  distance 

cl 

damage  function  dominates  the  circular  normal  distribution,  as 
illustrated  in  Figure  5A,B.  Notice  that  the  distance  damage  func- 
tion decreases  more  gradually  for  a larger  WR.  This  causes  the  cir- 
cular normal  distribution  to  be  centered  over  a smaller  and  smaller 
portion  of  the  damage  curve.  Consequently,  PD  goes  to  P,  as  WR/CEP 
increases  and  D/WR  remains  constant.  This  effect  is  evident  in  the 
tabular  results. 

3.4.2  Limit  of  PD  Along  Vertical  Axis 

Consider  the  limit  of  PD  as  the  offset  distance  D increases 
while  WR/CEP^  remain  constant.  Since  WR  remains  constant,  D/WR 
increases  as  the  offset  distance  increases.  Eventually,  the  aim- 
point  A becomes  so  far  offset  from  the  target  that  the  weapon  does 
not  interact  with  the  target  as  shown  in  Figure  5C.  In  other  words, 
the  probability  of  damage  to  the  target  is  0;  that  is,  PD  goes  to 
0 as  D/WR  increases  while  WR/CEP,  remains  constant.  Graphically, 
as  D increases  the  distributed  values  become  separated  so  the 
probability  of  damage  to  the  target  decreases  and  becomes  0. 
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